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Lattice kinetic equations incorporating the effects of external/internal force fields via a 
shift of the local fields in the local equilibria, are placed within the framework of contin- 
uum kinetic theory. The mathematical treatment reveals that, in order to be consistent 
with the correct thermo-hydrodynamical description, temperature must also be shifted, 
besides momentum. New perspectives for the formulation of thermo-hydrodynamic lat- 
tice kinetic models of non-ideal fluids are then envisaged. It is also shown that on the 
lattice, the definition of the macroscopic temperature requires the inclusion of new terms 
directly related to discrete effects. The theoretical treatment is tested against a controlled 
case with a non ideal equation of state. 



1. Introduction 

Lattice implementations of discrete-velocity kinetic models have gained considerable 
interest in the last decade, as efficient tools for the theoretical and computational investi- 
gation of the physics of complex flows (Shan & Chen (1993); Shan & Doolen (1995); Swift, 
Osborn & Yeomans (1995); Wagner & Pooley (2007); Li & Wagner (2007)). An impor- 
tant class of discrete- velocity models for ideal fluid flows, the lattice Boltzmann models, 
(Wolf Gladrow (2000); Benzi, Succi & Vcrgassola (1992); Chen & Doolen (1998)) can 
be derived from the continuum Boltzmann (BGK) equation (Bathnagar, Gross & Krook 
(1954)), upon expansion in Hermite velocity space of the single particle distribution func- 
tion, f(x, £, t), describing the probability to find a molecule at space-time location (x, t) 
and with velocity £ (Shan & He (1998); Martys, Shan & Chen (1998); Shan, Yuan & 
Chen (2006); Grad (1949)). As a result, the corresponding lattice dynamics acquires a 
more systematic justification in terms of an underlying continuum kinetic theory. 
The main goal of this paper is to extend such systematic link between continuum and 
lattice formulation to the case of hydrodynamical thermal fluctuations in non-ideal fluids 
under the action of both external (say, gravity) and/or internal forces. In particular, we 
shall show that if the effects of the force field are taken into account via a uniform shift of 
the momentum in the equilibrium distribution, as proposed in Shan & Chen (1993), the 
evolution for the total kinetic energy needs to be corrected as well. A viable possibility 
is to introduce an ad-hoc shift in the temperature field entering the local equilibrium. By 
doing so, the total kinetic energy recovers the correct thermo-hydrodynamical evolution. 
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Let us consider the usual continuum, D-dimensional, Boltzmann BGK equation 
%H- Vf+B*tf= ~ (/-/ (0) ) ; f^PAu) = J^e-X-^™ (1.1) 

where g represents an acceleration field and r a relaxation time towards the local equi- 
librium /(°). This local equilibrium will depend on the local density, p, momentum, pu, 
and temperature, 9. 

We will show that, as far as the macroscopic evolution of the hydrodynamical fields 
is concerned, it is possible to renormalize the action of the force pg only in terms 
of a suitable shift of the local Maxwellian equilibrium distribution appearing in (1.1): 
f( ^ (£; p, 9, u) — > /(°) (£; p, 9, u). The new -shifted- Boltzmann formulation being: 

| + * ■ V/ = ~{f - / (0) ); / (0) (€; P, 0, u) = ^-^e-l^l'/*. (1.2) 

The first result of the paper is to show that, in order to recover the correct thermo- 
hydrodynamic behaviour, the shifted local velocity and temperature must take the fol- 
lowing form: 

u = u + rg 9 = 9-r 2 g 2 /D. (1.3) 

The idea of shifting momentum has been pioneered for the case of a driving force due to 
internal, self-consistent interactions, which depend on the system configuration, typically 
the density field distribution. Such kind of extension has proven instrumental to the 
succesfull formulation of lattice kinetic theory of isothermal non-ideal fluids (Shan & 
Chen (1993); Shan & Doolen (1995)), in which one is concerned with the isothermal 
hydrodynamic evolution of density and momentum alone. 

Here we extend it further to the important case of thermal hydrodynamic fluctuations. 
Indeed, the use of shifted-equilibria has many immediate and important methodological 
consequences: (i) it provides an elegant way to incorporate the force effects thereby 
dispensing with the need of taking derivatives of the distribution function in velocity 
space; (ii) it allows a systematic derivation of lattice kinetic equations for non-ideal 
fluids with (pseudo)-potcntial energy interactions (Shan, Yuan & Chen (2006); Martys, 
Shan & Chen (1998); Guo, Zheng & Shi (2002)); (iii) it highlights the need of including 
a suitable redefinition of the hydrodynamical fields on the lattice, in order to recover the 
correct continuum limit of the thermo-hydrodynamical equations for density, momentum 
and total kinetic energy, K. Here and throughout, by continuum thermo-hydrodynamic 
limit we refer to the following set of macroscopic equations (repeated indices are summed 
upon): 

d t p + di{pui) = 

dt(pu k ) + di(Pi k ) = pgk (1-4) 
d t K + \diqi = pg t Ui 

where Pik and qi are momentum and energy fluxes (still unknown at this level of de- 
scription). More precisely, we shall show that the above equations can be obtained ex- 
actly from the previous continuum Boltzmann Equation with shifted equilibrium (1.2). 
This would not be the case, if the momentum is the only shifted quantity in (1.2). The 
relevance of such a continuum kinetic theory is mainly motivated by the final goal of 
formulating lattice versions of the Boltzmann Equations for non-ideal fluids, including 
thermo-hydrodynamic effects. In fact, in lattice formulations, the need of representing 
velocity degrees of freedom through a limited set of discrete speeds, raises the problem 
of a correct and efficient implementation of the continuum velocity-streaming operator 



Lattice Boltzmann method with self- consistent equilibria 



3 



g ■ V^/. We shall show that shifted equilibria in the continuum representation have a 
"well-defined lattice analogue, so that the lattice counterpart of the continuum description 
(1.2) can be obtained through the usual lattice Boltzmann discretization: 

fi(x + c,At,t + At)-f l (x,t) = - — (Mx,t)-f\x,p^M L \e^)) (1.5) 

T 

where the subscript I runs over the discrete set of velocity on the lattice, C;, and the 
superscript L indicates that the macroscopic fields are now defined in terms of the lattice 
Boltzmann populations: 

p {l) = E /<; p {L)uiL) = E D p {L)diL) = E ^ u(L) \ 2 fi- 
i i i 

In the expression (1.5), the Boltzmann equilibrium (sec Nie, Shan & Chen (2008) for its 
explicit expression) is computed with shifted momentum and temperature, as follows: 

After some algebra, it can be shown that the temperature shift AO can be expressed in 
closed form as a function of the lattice time-step At: 

Moreover, in order to recover the thermo-hydrodynamical equations 

'd tP W +d i {puf ) )=0 

< dt(p (H) u[ H) ) + d t (Pi H} ) = pWg k (1.6) 

the hydrodynamical fields can be computed in terms of a closed expansion to all orders in 
At and can be calculated in terms of a suitable lattice operator. For example, density is 
left unchanged, 

P {H) = P, while the first non trivial correction is given by the well-known 
pre and post-collisional momentum average (Buick & Greated (2000)): 

U {H) = U {L) A* 

2 y 

as well as by a new, non-trivial, correction to the temperature field: 

(H) (L) (At)V 
4D ' 

First, we notice that in the limit A< — > 0, the lattice formulation for both shifted fields 
and hydrodynamical fields goes back to the continuum one, as it should. Second, the 
continuum formulation in terms of shifted fields indicates a straightforward link with the 
discrete variables via the Hermite-Gauss expansion (Shan, Yuan & Chen (2006)). Third, 
and maybe more important for applications, we emphasize that in order to achieve a 
self-consistent thermo-hydrodynamical description in the lattice (1.6), both momentum 
-as it was known- and temperature, acquire discrete corrections. 

The paper is organized as follows: in section 2 we present the basic ingredients of the 
continuum BGK model with shifted equilibria and we show how to compute the correct 
normalization of the local fields, so as to recover the correct macroscopic equations. 
This procedure is extended to the lattice models in section 3 with a detailed analysis 
of the discrete contribution to the rcnormalization procedure. A numerical test of the 
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above arguments is provided in section 4. Finally, in section 5, we will briefly discuss 
the physical interpretation of the shifted continuum model and propose some further 
development. 



2. Shifted Continuum Equilibrium 

In this section wc deal with the macroscopic properties of a continuum model (1.2). 
The main goal is to renormalize the effects of the forcing term g ■ V^/ in (1.1) via a 
suitable local equilibrium with shifted fields: 

u = u + Au(g,T); = + A%,r). (2.1) 

It is well known (Bathnagar, Gross & Krook (1954); Wolf Gladrow (2000)) that the usual 
definition of Boltzmann Equation with explicit forcing given in (1.1) leads to the exact 
macroscopic equations (1.4) with the averaged fields given by 

P = J «/; pu = J d^f; K. = 1 -{ P D9 + pu 2 ) = lj dtt 2 f- (2.2) 
The momentum and energy fluxes 

Pij = J debt,/ q t = J dt^M 

in the LHS of equations (1.4) need to be closed, a task which is usually accomplished via 
the Chapman-Enskog expansion. In order to derive (1.4) from (1.1) it is indeed sufficient 
to notice that the collision operator —\{f — f^) preserves both mass, momentum and to- 
tal kinetic energy, as long as the local equilibrium is expressed in term of the macroscopic 
fields, p,u,6, i.e. whenever, besides the definition (2.2), the local equilibrium carries the 
same mass, momentum and energy of the actual distribution, namely: 

J d£/ (0) ; pu = J d^/ (0) ; K. = \{ P D6 + pu 2 ) = \J ^ 2 / (0) - (2.3) 

When considering the effect of the forcing field embedded within the shifted equilibrium 
(1.2), the only difference with the previous standard derivation comes from the fact that 
now the averaged macroscopic fields, when evaluated on the shifted equilibrium, do not 
coincide with the hydrodynamical fields defined in terms of the local particle distribution 
(2.2), i.e. the collision operator, — -(/ — f^) preserves momentum and energy only 
globally, but not locally. Still, it is easy to realize that the extra momentum and energy 
brought by the shifted equilibrium is given by: 

pu-pu = J dtt(I {0) -f); \(pDe+pu 2 )- l -{ P D6+pu 2 ) = l -J dee 2 (/ (0) -/) (2.4) 

and that if we chose the shifted fields as given by expressions (1.3) the exact macroscopic 
equations (1.4), for density, momentum and kinetic energy evolution, are recovered. To 
this purpose, it is sufficient to evaluate the first three lowest momenta of eq. (1.2) and 
use the relations (2.4) and (2.2). This is the first result of the present work. 
Let us stress once again that the shift in the temperature only responds to the need of 
cancelling out extra terms in the rhs of eq. (1.4), that would otherwise result from the 
momentum-shift alone. One may wonder if beside the formal correct unclosed equations 
(1.4) the two Boltzmann Equations formulation (1.1) and (1.2) do also share the same 
hydrodynamical behaviour, i.e. if the unclosed momentum and heat fluxes do have the 
same closure. By performing the whole Chapman-Enskog expansion, it can be shown that 
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this is indeed the case, at least up to second order in the expansion parameter where 
dissipative terms in momentum and heat appear. Details of these calculations will be 
reported elsewhere. 



3. Lattice Implementation 

In this section, we treat the lattice averaged equations and discuss the way that the 
corresponding shifts in the momentum and temperature fields are affected by the dis- 
cretization of the algorithm. To this purpose, let us go back to the lattice Boltzmann 
equation with shifted momentum and temperature given by expression (1.5). First, we 
Taylor expand the LHS of (1.5), and obtain the lattice-series expression 

with Di t = dt + ci ■ V. We can then rewrite (1.5) in a compact form 

(e AWw - 1) fi = Atd (3.1) 

with Ci = —- (fi — fi J the collisional operator. A formal inversion of (3.1) yields: 



_ AtD ltt _ 



AtD Lt e- AtD >>* i 



(! _ e -AtA.t) 



Ci 



where the LHS is recognized as the generating function of Bernoulli polynomials B n (x) 
(Gradshtein & Ryzhik (2000)) : 



n-l 



n=0 



Let us also introduce the operator A, At, defined by the Taylor expansion in DuAt: 



DitCi.At = 



AtD u e 



-At A 



(1-e 



-AtD, 



- 1 



AtD n { AtDn (AtDit) 3 {AtDitf ,, A 

' 1 ^ + ^ ^- - ^ ^- + O { (AtD, 

6 360 15120 vv ' 



The relevant point is that the above operator can be rewritten in terms of a lattice 
operator performing an inverse shift in space and time: S^a* = (e _A * £>! t — 1). It is easy 
to realize that its action on any field, say </>, defined on the lattice gives back: 

Si A t<t>{x, t) = 4>(x - Atc u t - At). 

The action of the operator C^At can be recast as follows: 

A. At = I 1 + ^Sl.At — T^S? A f + — — Sf A t — —Si At. + O (S? 



LAt - — ^ + g^At ~ ^LAt + ^,At ~ ^LAt + ^ ^,At, 

This shows that it is possible to rewrite the final dynamics, exactly to all orders in At, 
by retaining only shift-operators on the lattice: 

A,t {f i - £i,AtCi) = Ci 

Alternatively, by recalling the definition of C/ , we also have: 



D 



Li 



( K fl + ^ l ,A t (f l -f[° ) ))=- 1 -(f l -fh- 



We can then average this equation in velocity space and look at the equations for the first 
three-order momenta, so as to recover the hydrodynamic evolution for density, momen- 
tum and total energy (1.6). Simple calculations shows that the LHS set of macroscopic 
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equations (1.6) is obtained by means of the following definitions of macroscopic hydro- 
dynamic (H) fields (scalars, vectors and second-order tensors) 

= 4ft + i (e, cjCtMfi - // 0) 

— J, J t i i/v J J r, . At, _ 



= (% P ew + \p{u^f) = i e, <f/« + i (i <f A,A t [/* - /, (0) ]) 



af 10 = I E, c?cft + H i E ; c?4a,a* [/i - /T: 



In order to capture the correct RHS of (1.6) as well, we must choose the fields entering 
the shifted equilibrium, /( 0) in (1.5), = m< l) +Am, 8^ = 6»( L '+A6>, in the following 
form: 

{Am = rg 
A6 = %(pg i {u ( f ) -u i )-^f). 

Let us notice that the above expression for the temperature shift is implicit, i.e. it is given 
in terms of the hydrodynamical velocity vS H > which depends itself on the equilibrium. 
One may get a closed expression only via a Taylor expansion in At. For example, to first 
order in the expansion of £i,At, we have simply 

corresponding to the following temperature shift: 

= 9 (L) + r(At- r)g* 

Consequently, the hydrodynamical velocity and temperature become: 

f = + fg 

where the hydrodynamic velocity is nothing but the pre and post-collisional average, 
while the non trivial correction in At to the temperature fields is the new result, as 
anticipated in the introduction. 



4. Numerical test 

We now proceed to the implementation of a test case of the above procedure, where the 
need for the temperature shift appears in full. The most important instance where total 
energy conservation is crucial is the case of a gas(fluid) departing from ideal conditions, as 
a result of an internal, self-consistent, potential. In this case, the thermo-hydrodynamical 
equations must conserve the total energy, given by the sum of the total kinetic plus the 
potential energy. Typical relaxation experiments will then show a non trivial exchange 
between the kinetic and potential energy components, until a dynamical or static equi- 
librium is finally attained. 

We specialize the discussion to a simple, and yet non-trivial, case where the interpar- 
ticlc force is purely repulsive : 

F = pg(x,t)=g P (x,t)J2wip(x + c l At,t)ci « -VP h (mt) P^ip) = ~\0p 2 (4.1) 

i 



with wi a suitable set of weights which enforce the right symmetries on the lattice (Shan 
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& Chen (1993); Shan & Doolcn (1995)). Repulsion is imposed by choosing a negative 
coupling constant, Q < 0. This case allows full control of the non-ideal part of the 
equation of state. Indeed, the bulk pressure provides the usual ideal-gas contribution, 
■f& — p0 plus the non-ideal one, given by the Taylor expansion of the forcing term in 
(4.1), p^ tot ^ = p( id > -\-p( tnt '. The system also has internal potential energy Ey — ~^Gp 2 , 
where we have neglected possible contributions coming from strong density gradients. The 
transport equation for this intermolecular potential energy (He & Doolcn (2001); Snider 
(1991, 1995)) reads as follows: 

d t E v + d k {u[ H) E v ) = -(a^fVf", 

as can be readily derived from the density evolution in (1.4). By summing to the total 
kinetic energy evolution, we obtain the following total energy balance: 

d t (Ev + = -^(.f'pf) - \d4 H) - d k (ui H) Ey), 

whose divergence-like structure at the RHS, ensures total energy conservation. 

Let us stress that the reconstruction of a total divergence is only possible thanks to 

the superposition of the contribution (djUj')P^ , stemming from the evolution of the 

potential energy, plus the contribution Uj H \dj P^ mt ^ ) , coming from the RHS of the total 
kinetic energy in (1.6). Here, we fully appreciate the importance of the temperature shift, 
in order to recover the correct total energy dynamics. Using fully periodic boundary 
conditions, the shifted lattice Boltzmann formulation is therefore expected to provide 
conservation of the total energy from the hydrodynamical point of view. 

To ensure a sufficiently accurate recovery of the thermal transport phenomenon, we 
employ a two-dimensional 37-speed LB model, corresponding to a ninth-order accurate 
Gauss-Hermitc quadrature. In conjunction, the following fourth-order Hermitc expansion 
of the Maxwellian is used as an equilibrium distribution (Shan, Yuan & Chen (2006); 
Nie, Shan & Chen (2008)). A simulation is performed on a L x x L y = 10 x 100 grid, with 
a small perturbation of a single sinusoidal wave in the temperature field 0( H \x,y,t = 
0) = 1.0 + esin(27rj//L s/ ) (e = 0.01). The initial density field is constant. The difference of 
the total energy J ((Ey + K^ H ^)){t)dxdy with respect to its initial value is monitored in 
figure 1. This figure clearly shows that the lattice Boltzmann formulation without shifted 
temperature is not able to sustain satisfactory energy conservation. On the other hand, 
upon shifting the temperature field, the correct energy balance is recovered. The energy 
conservation is still below machine precision, due to the fact that our expression for the 
total energy is given in terms of a continuum description of the non-ideal forcing term 
(4.1). In order to further improve the accuracy of energy conservation, a discrete version 
of the internal potential energy, Ey , thermodynamically consistent on the lattice, needs 
to be developed. 



5. Conclusions and Perspectives 

Before concluding, let us discuss further the physical meaning of the shifted Boltzmann 
Equation (1.2). It is easy to realize that it can also be rewritten as: 

df{X dt tt] + £ • Vf(x, €,«) = — (/(*, £, t) - e -ff-v t -T a £v t .v t/ (o) (€ . pj 9j u) \ , 
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Figure 1. Variations of total energy for a non ideal system with an initial sinusoidal wave 
in the temperature field 6 t ~ H \x,y,t = 0) = 1.0 + e sm{2iry/L v ) with e = 0.01 and smooth 
hydrodynamical velocity fields. The lattice Boltzmann parameters in (4.1) and (1.5) are such 
that r/ At = 0.6, Q — —3.0. Two simulations are carried out. The first simulation is only with 
shifted momentum (A). A second one with shifted momentum and temperature (o). 



Upon Taylor expanding the RHS terms corresponding to the shift in the momentum. 
— e~ T9 ' v «_f(°) . up to second order in r, we obtain: 

% + € • V/ = -± (/ - /C°>) - g • V c /<°> + T -gg : V £ V £ /<°> (5.1) 

It is simple to check that, as for the isothermal dynamics of density and momentum, the 
shifted equilibrium Boltzmann equation (1.2) is equivalent to (5.1). A stabilizing diffu- 
sion term in velocity space stands therefore out. 

It is also noted that, as a first order of approximation, this kind of diffusion term can be 
thought as deriving from a standard BGK dynamics, with the equilibrium distribution 
f(°\x,£,t) replaced by a smoothed version, resulting from coarse-graining in velocity 
space, filtering fluctuations up to 5v' < gr. Considering the Taylor expansion of the 

temperature shift operator, — e _T iD V «' v « f(°\ we obtain, up to r 2 , an extra-term pro- 
portional to: 

5 2 A,/(°). (5.2) 

Such contribution can be interpreted as deriving from a stochastic component in the 
acceleration field. To illustrate the point, let us start again from the continuum BGK 
equation (1.1) and let us consider the streaming term in velocity space, -4 • Vf /, in which 
the molecular velocity & obeys the following Langcvin equation: 

H = m 

where r\i is a standard delta-correlated Gaussian noise with zero average and the following 
normalization: 

v7i) = (rurjj) = ^aSij. 

We expect the stochastic term to provide a mechanism for producing thermal fluctuations 
in the fluid. The key ingredient is a correct evaluation of the term (^V^/) where £i 
is now a stochastic variable and (...) stands for an average over possible realizations of 
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the stochastic term. We can apply Novikov's theorem (Novikov (1964)) and approximate 
& ~ i!t T = Vi T - This yields: 

(^VeJ) = (^6-)(v ft v^/> « T-taXVfrVe,/) 

which is of the same form of (5.2) for a suitable choice of cr. Let us emphasize that 
the shifted equilibrium formulation (1.2) is equivalent to the standard BGK formulation 
(1.1), as far as the macroscopic equations (1.4), and their Chapman-Enskog expansion 
up to second order, are concerned. The diffusive extra-terms stemming from the Tay- 
lor expansion indicate that, although maintaining the same hydrodynamics, the shifted 
equilibrium formulation should nonetheless feature better stability with concern to the 
global evolution of the probability density, f(x,£,t). 

Summarizing, we have investigated lattice kinetic equations incorporating the effects 
of extcrnal/sclf-consistent force fields via a shift of the local fields in the local equilibria. 
The mathematical treatment reveals that, besides momentum, temperature also receives 
a self-consistent shift from the force field. The contribution of the temperature shift can 
also be traced back to a stochastic component in the acceleration field, thus pointing 
to potentially new directions for the formulation of lattice Boltzmann models for non- 
ideal fluids with thermo-hydrodynamic transport effects (Brennen (2005); Rowlinson & 
Widom (1982); De Gennes, Brochard-Wyart & Quere (2003)). Work along these lines is 
in progress, which will hopefully permit to attack a broad class of complex flow problems 
with thermal effects, such as thermally driven phase-transitions, crystal growth, melting 
and many other non-equilibrium thermo-hydrodynamic transport problems. 

Finally, it has been shown that, in order to recover the correct thermo-hydrodynamical 
equations (Boghosian (2008)) on the lattice, the macroscopic temperature must acquire 
new terms, directly related to the lattice spacing. These terms naturally vanish in the con- 
tinuum limit, thus preserving the consistency of the discrete theory. Many items remain 
open to future investigation. For instance, it would be interesting to extend the present 
treatment to more general collision operators, including multi-time relaxation models 
(Dellar (2003)), which would permit to model fluids at non-unitary Prandtl numbers 
(Ahlers, Grossmann, & Lohse (2009)). The establishment of a H-theorem for continuum 
and discrete kinetic equations with self-consistently shifted equilibria, also appears an 
interesting topic for future research (Succi, Karlin & Chen (2002); Ansumali & Karlin 
(2005)) 

REFERENCES 

Ahlers C, Grossmann S., and Lohse D. 2009 Heat transfer and large-scale dynamics in 

turbulent Rayleigh-Benard convection. Rev Mod. Phys. in press. 
Ansumali S & Karlin I. 2005 Consistent Lattice Boltzmann Method. Phys Rev. Lett., 95, 

260605. 

Bathnagar P.-L., GROSS E. & Krook M. 1954 A model for collision processes in gases. 

Physical review 94, 511-525 
Benzi R., Succi S. & Vergassola M. 1992 The lattice Boltzmann equation: theory and 

applications. Phys. Rep. 222, 145-197 
Boghosian B. M. 2008 Exact Hydrodynamics of the Lattice BGK Equation. arXiv:0810.2344vl 
Brennen C. 2005 Fundamentals of Multiphase Flow. Cambridge University Press 
BuiCK J. M. & Created C. A. 2000 Gravity in a lattice Boltzmann model. Phys. Rev E 61, 

5307-5320 

Chen S. & Doolen G. 1998 Lattice Boltzmann method for fluid flows. Annu. Rev. Fluid Mech. 
30, 329-364 



10 M. Sbragaglia, R. Benzi, L. Biferale, H. Chen, X. Shan and S. Succi 

De Gennes P. -J., Brochard-Wyart F. & Quere D. 2003 Capillarity and Wetting Phenom- 
ena. Springer, New York 

Dellar P. J. 2003 Incompressible limits of lattice Boltzmann equations using multiple relax- 
ation times. Jour. Comp. Phys. 190, 351-370 

Grad H. 1949 Note on N-dimensional Hermite polynomials. Pure Appl. Math. 2, 325 

Gradshtein I.S. & Ryzhik I.M. 2000 Tables of integrals and series. 6th edition, Academic 

Guo Z., Zheng C. & Shi B. 2002 Discrete lattice effects on the forcing term in the lattice 
Boltzmann method, Phys. Rev. E 65, 046308 

He X. & Doolen G. 2001 Thermodynamic Foundations of Kinetic Theory and Lattice Boltz- 
mann Models for Multiphase Flows. Jour. Stat. Physics 107, 309-328 

Li Q. & Wagner A. J. 2007 Symmetric free-energy-based multicomponent lattice Boltzmann 
method. Physical Review E 76, 036701 

Martys N. S., Shan X. & Chen H. 1998 Evaluation of the extrenal force term in the discrete 
Boltzmann equation. Phys. Rev. E 58, 6865 

Nie X., Shan X. & Chen H. 2008 Thermal lattice Boltzmann model for gases with internal 
degrees of freedom. Phys. Rev. E 77, 035701(R) 

Novikov E. A. 1964 Functionals and the method of random forces in turbulence theory. Zh. 
Exper. Teor. Fiz. 47, 1919-1926 

Rowlinson J. R.& Widom B. 1982 Molecular theory of Capillarity. Clarendon, Oxford, 1982 

Shan X. & Chen H. 1993 Lattice Boltzmann model for simulating flows with multiple phases 
and components. Phys. Rev E 47, 1815 

Shan X. & Doolen G. 1995 Multicomponent lattice-Boltzmann model with interparticle in- 
teraction. Jour. Stat. Physics 81, 379-393 

Shan X. & He X. 1998 Discretization of the Velocity Space in the Solution of the Boltzmann 
Equation. Phys. Rev. Lett. 80, 65-68 

Shan X., Yuan F. & Chen H. 2006 Kinetic theory representation of hydrodynamics: a way 
beyond the NavierStokes equation. Jour. Fluid Mech. 550, 413-441 

Snider R. F. 1991 A reinterpretation of dense gas kinetic theory. Jour. Stat. Phys. 63, 707-718. 

Snider R. F. 1995 Conversion between kinetic energy and potential energy in the classical non 
local Boltzmann equation. Jour. Stat. Phys. 80, 1085-1117. 

Swift M. R., Osborn W. R. & Yeomans J. M. 1995 Lattice Boltzmann Simulation of 
Nonideal Fluids. Phys. Rev. Lett. 75, 830-833 

Wagner A. J. & Pooley C. M. 2007 Interface width and bulk stability: Requirements for the 
simulation of deeply quenched liquid-gas systems. Physical Review E 76, 045702(R) 

Wolf-Gladrow D. 2000 Lattice-Gas Cellular Automata And Lattice Boltzmann Models. 
Springer, New York, 

Succi S., Karlin I. & Chen H. 2002 Colloquium: Role of the H theorem in lattice Boltzmann 
hydrodynamic simulations. Review of Modern Physics Colloquia 74, 1203 



